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I.  INTRODUCTION 


When  a  self-pinched  particle  beam  propagates  through  a  gas.  It 
undergoes  multiple  scattering  off  background  nuclei.  As  a  result,  the  mean 
radius  of  the  beam  tends  to  Increase  and  the  mean  beam  density  decreases. 
Recent  experiments  at  Los  Alamos  National  Laboratory  (LANL)  using  the 
PHERMEX  electron  beam  have  called  into  question  standard  formulae  used  to 
compute  the  rate  of  expansion  of  the  beam.  In  particular,  the  length  for 
one  e-folding  in  mean  radius,  the  so-called  Nordsieck  length,  was  measured 
to  be  about  a  factor  of  two  longer  than  the  theoretical  prediction.  In 
this  report,  we  show  that  this  discrepancy  is  due  to  inaccurate  treatment 
of  small  angle  multiple  scattering.  Section  II  describes  a  more  accurate 
treatment  based  on  the  Moliere  formulation  of  multiple  scattering.  We  show 
how  a  numerical  algorithm  was  constructed  to  implement  this  improved  treat¬ 
ment  in  a  particle  simulation  code.  In  Sec.  Ill,  we  give  the  results  of 
numerical  computations  of  Nordsieck  length  as  a  function  of  beam  power, 
using  a  simple  particle  code  SCATR.  On  the  basis  of  these  results,  we 
carry  out  a  semi-analytical  derivation  of  an  expression  for  the  Nordsieck 
length.  Section  IV  compares  numerical  and  analytical  results  with  experi¬ 
mental  data  for  the  PHERMEX  and  ASTRON  electron  beams.  Good  agreement  is 
obtained.  In  Sec.  V,  we  examine  the  evolution  of  a  non-Bennett  equilibrium 
(an  annular  beam)  and  show  that  it  becomes  Bennett-like  in  less  than  one 
Nordsieck  length.  We  summarize  our  conclusions  in  Sec.  VI. 
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II.  THEORY  OF  MULTIPLE  SMALL  ANGLE  SCATTERING 


When  a  charged  particle  passes  through  a  medium  made  up  of 
scattering  centers  with  cross  section  <j(n)d£i,  the  number  of  scattering 
events  which  it  undergoes  in  traveling  a  distance  I  is 

N0  =  n tf  o(n)dn 
00 

=  2nnt  /  o(e)ede  ,  (1) 

o 

where  n  is  the  number  density  of  scattering  centers,  and  e  is  the  angular 
deflection.  In  writing  Eq.  (1)  we  have  made  the  assumption  that  the  net 
deflection  of  the  incident  particle  is  small  (much  less  than  a  radian)  so 
that  the  path  of  the  incident  particle  is  approximately  a  straight  line. 
We  have  also  assumed  that  the  cross  section  is  azimuthally  symmetric.  When 
N0  >  1,  we  are  in  the  multiple  scattering  regime. 

A.  Fokker-Planck  Formalism 

To  describe  the  behavior  of  a  large  number  of  incident  particles, 
we  use  a  distribution  function  f(r^,p^,t)  where  ,  p^  are  the  transverse 
phase  space  coordinates  of  the  beam  of  particles.  The  longitudinal  motion 
of  the  particles  is  assumed  to  be  given  by  z  =  vt  where  v  •  c  is  the  bean 
velocity.  The  distribution  function  obeys  the  transport  equation 


=  -2*nv  f  /  o(x)xdX 

+  2*nv  /  f(e  -  x)  o(x)  dx  (2) 


2 


dp, 

where  is  the  force  due  to  self  fields,  e  is  a  vector  in  the  transverse 
plane  related  to  p^  by  e  =  P^/ymc,  m  is  the  electron  mass  and 
y  •  (l-v2/c2)-1/2.  For  many  purposes,  the  Fokker-Pl ancle  treatment  of  the 
collision  operator  on  the  right-hand  side  is  adequate.  This  consists  of 
an  expansion  to  second  order  in  x  yielding 


/af  \  .1 

/*2\ 

2  "I 

3  t  l  if 

Wcoll  '  T 

<*  > 

362  6  36 

where 


2irnv 


3 

/  x  dx  o(x) 
0 


(3) 


The  small  angle  Rutherford  cross  section  for  an  electron  scattering  off  a 
nucleus  is 


2iro  (x)  xdx 


8irZ(Z+l) 

A  2 
e  y 


x 


2“  Xdx/X 


(4) 


where  -e  is  the  electron  charge,  Ze  the  nuclear  charge,  and  &  =  v/c.  Sub¬ 
stituting  this  into  Eq.  (3),  we  find  that  the  integral  for  <x2>  is  improper 
at  both  the  lower  and  upper  limits  of  Integration.  The  divergence  at  the 
lower  limit  is  due  to  the  long  range  of  the  bare  Coulomb  potential.  A 
suitable  cutoff  can  be  introduced  by  taking  account  of  the  screening  effect 
of  atomic  electrons.  For  the  Thomas-Fermi  atomic  model,  for  example,  one 
obtains1  the  lower  bound  6min  =  */a  =  */(0.855aoZ-1/3)  where  *  is  the 
deBroglie  wavelength  of  the  incident  electron,  and  a0  is  the  Bohr  atomic 
radius. 
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The  divergence  at  the  upper  limit  is  an  artifact  of  the  small 
angle  approximation  used,  since  a  particle  cannot  scatter  through  an  angle 
greater  than  it.  The  Fokker-Planck  equation  can  be  derived  only  in  the 
small  angle  limit,  however,  so  that  one  can  either  (a)  introduce  some  reas¬ 
onable  upper  bound  on  the  scattering  angle  or  (b)  go  to  a  different  treat¬ 
ment  of  small  angle  scattering.  Two  commonly  used  upper  bounds  are 
emax  a  n  f°r  l°w  electron  energies  (y  <  100)  and  e^ax  =  */Rn,  where 
Rn  is  the  nuclear  radius,  used  by  Rossi  and  Greisen,  for  high  energies2 
(y  >  100).  Using  either  of  these,  one  can  solve  Eq.  (2)  analytically  for 
certain  initial  conditions.1*3  For  future  reference,  we  give  the  results 
for  two  cases: 

(i)  Noninteracting  Particles 

Assume  the  incident  beam  to  have  infinite  radial  extent  with 
3/arj^  =  0  and  let  f  «  6(v^)at  t  =  0.  The  solution  to  Eqs.  (2)  and  (3)  is 

then 

f(e,0  -  exp[-e2/<x2>£]  ,  (5) 

where  i  =  vt.  Since  <i2>  =  2na2  An(emax/0min)  is  a  constant,  the  beam 
acquires  a  Maxwellian  distribution  in  transverse  velocity  whose  temperature 
increases  linearly  with  the  distance  propagated. 

( i i )  Self-Pinched  Bennett-Profi  •>  Beam 

From  Eqs.  (2)  and  (3),  Lee4  has  shown  that  a  Bennett  equilibrium 
beam  expands  self-similarly  due  to  scattering  with  Bennett  radius 

R(£)  =  Rq  exp(£/LN)  (6) 
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where  is  the  so-called  Nordsieck  length.  Equation  (6)  is  valid  as 
long  as  v/Lr  «  where  wg  is  the  mean  betatron  frequency,  i.e., 

the  beam  expansion  is  quasi -static.  The  expression  obtained  for  is 

ln  =  ;  2 '2  *n(em;7e~T  *  (7) 

ny  a  '  max  min' 

where  v  is  Budker's  parameter.  Here  and  subsequently,  v  is  a  measure  of 
the  net  current  flowing.  The  plasma  current  is  assumed  to  have  the  same 
profile  as  the  beam  current.  Taking  emax  =  k/R^.  we  obtain  £n(emax/ 
©min)  "  2zn(210Z-1/ 3 ) .  For  0fnax  =  7T»  ^n( ©max/ ©min)  =  An(192yeZ  "*/") 
We  shall  see  that  these  values  for  emax  overestimate  the  effect  of  small 
angle  scattering  over  a  wide  range  of  parameters  relevant  to  present-day 
beam  propagation  experiments. 

B.  Moliere  Formalism 

For  the  case  of  noninteracting  particles  with  df/dr^  =  0  and 
f (t=0)  «  6 ( ) ,  Eq.  (2)  can  be  solved  exactly  for  arbitrary  o(e)  by  the  use 
of  Fourier-Bessel  transforms.  This  was  first  done  by  Moliere5,  but  the 
brief  description  given  here  is  from  a  simpler  derivation  of  Moliere' s 
equation  by  Bethe.1  The  solution  takes  the  form 

f(e,t)  =  /  n  dn  J  (ne)  eN^  '  No  ,  (8a) 

o 

00 

where  N(n)  =  n i  J  o(x)JQ(nx)  X  dX>  (8b) 

o 

and  J0  is  the  Bessel  function  of  order  zero. 
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When  N0  >  1,  only  values  of  n  for  which  N0  -  N{ n)  <  N0  contribute  in 
Exj .  (8a).  For  this  regime,  it  turns  out  that  Eq.  (8b)  can  be  evaluated 
analytically.  To  do  this  most  conveniently,  one  changes  variable  to  c  = 
n(nto2B)1/2,  where  B  is  the  solution  to 

eB/B  =  n*a2/(1.32e^in)  -  yi.32  .  (9) 

Then  one  can  write 

2  2  2 

n(o  -  n0  -  -  f-  +  %  *n  r  (10) 

for  0  <  c  <  B1/2.  Substituting  this  into  Eq.  (8a),  one  obtains  Moliere's 
expansion  in  1/B: 

f (e)ede  -  ?dc[f(0)U)  +  B'1  f(1)U)  +  B"2  f(2)U)  +  ...]  (U) 

where  ;  -  e/(n*a2B) x/2 ,  and  f(°)  =  exp(-c2).  Thus,  the  leading  term 
in  the  Moliere  expansion  is  a  Gaussian  of  width  6m  =  (ma2B)1/2.  When 
B  is  large  (>  10)  we  can  write  B  «  in  N0  <*  fcn(jt).  So,  in  contrast  to  the 
result  in  Eq.  (5),  we  find  that  the  transverse  temperature  of  the  beam 
increases  as  l  *n(fc)  instead  of  as  l. 

Comparing  the  widths  of  the  Gaussians  from  Eqs.  (5)  and  (11),  we 
find  that  Eq.  (5)  overestimates  the  mean  square  transverse  momentum  of  a 
beam  of  noninteracting  particles  propagating  in  air  (see  Fig.  1).  Because 
the  width  of  the  Moliere  Gaussian  increases  faster  than  linearly  with  i,  it 
eventually  catches  up  to  the  Rossi -Grei sen  value,  but  only  after  -  104 
meter-atmospheres. 

Having  derived  the  Moliere  result,  one  can  go  back  and  deduce 
what  value  one  should  use  for  emax  to  obtain  the  Moliere  Gaussian  width. 
One  finds  that  emax  =  9m  is  the  appropriate  cutoff,  i.e., 

eM 

X2  l  =  e^j  =  2irnt  /  x3dx  o(x) 

®min 
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This  leads  to  a  transcendental  equation  for  which  is  equivalent  to 
Eq.  (9).  The  main  effect  of  the  higher  order  terms  in  Eq.  (11)  is  to  pro¬ 
duce  a  long  tail  on  the  distribution.  For  large  values  of  c»  the  expres¬ 
sion  for  f  goes  over  to  the  Rutherford  single  scattering  form.  As  a  result 
of  this  tail  on  the  distribution  function,  one  finds  that  the  shape  of  the 
distribution  function  changes  as  a  function  of  propagation  distance,  as 
shown  in  Fig.  2.  The  greater  the  propagation  distance,  the  closer  the  dis¬ 
tribution  approaches  the  Rossi -Greisen  Gaussian,  and  the  smaller  the  per¬ 
centage  of  particles  in  the  tail. 

When  self-field  effects  are  retained,  Eq.  (2)  becomes  rather 
intractable.  Analysis  along  the  lines  of  Ref.  4  is  made  difficult  by  the 
fact  that  the  rate  of  transfer  of  energy  to  the  transverse  direction  is  not 
constant,  but  depends  on  the  scattered  state  of  the  beam.  In  Sec.  Ill  we 
show,  however,  that  a  semi-analytical  treatment  can  be  performed  for  the 
special  case  of  a  Bennett  equilibrium. 

C.  Numerical  Algorithm 


In  order  to  look  at  problems  where  self-fields  are  important,  we 
devised  a  numerical  algorithm  to  implement  Moliere  scattering  in  particle 
simulation  codes.  In  the  algorithm,  scattering  angles  are  picked  randomly 
from  a  Moliere  distribution.  The  distribution  used  is  one  corresponding  to 
a  distance  of  propagation  Li  much  less  than  the  characteristic  length  over 
which  the  beam  radius  changes  due  to  scattering.  To  compute  the  Moliere 
distribution  one  has  to  perform  time-consuming  numerical  integrations,  so 
we  divide  the  distribution  into  two  parts, 


c  <  c. 


P(c)  dc 


1  ili 
B  4 
C 


c  > 


1 


where  is  an  angle  at  which  the  Rutherford  expression  becomes  a 

good  approximation  to  the  exact  Moliere  expression.  At  ccut  =  6,  for 
example,  the  discrepancy  between  the  two  expressions  is  15%.  For  ?  < 
CC(Jt»  an  array  of  scattering  angles  is  computed  once  and  for  all  on  the 
first  pass  through  the  algorithm.  For  c  >  ?cut,  scattering  angles  can 
be  obtained  directly  from  an  analytic  expression,  and  it  is  thus  not  too 
time-consuming  to  generate  new  angles  each  time  a  scattering  is  required, 
especially  since  large  angle  scatterings  are  relatively  rare.  The  actual 
choice  of  a  scattering  angle  proceeds  as  follows.  A  random  number  X  is 
chosen  between  0  and  A  where  A  =  /  P(c)d?.  If  X  <  A^,  where  A^  =  / 
P(c)dc  then  the  angle  corresponding  to  X  is  obtained  from  the  permanent 
array.  If  X  >  A^,  the  corresponding  angle  is  computed  analytically. 

The  length  At  cannot  be  made  arbitrarily  small,  because  Moliere 
theory  breaks  down  rapidly  for  N0(At)  <  20.  To  avoid  some  numerical 
integration  problems,  we  generally  choose  At  such  that  N0  40  in  our 
algorithm.  When  a  beam  expands  to  the  point  where  -  v/L^,  this 
value  of  At  can  become  comparable  to  the  length  scale  over  which  the  beam 
radius  changes  (see  Sec.  IVB).  We  then  resort  to  the  treatment  of  Keil, 
Zeitler,  and  Zinn6»7  which  is  specifically  addressed  to  the  regime  1  j< 
N0  _<  20,  the  so-called  “plural"  scattering  regime.  We  can  thus  chose  a 
suitably  small  At  while  maintaining  accuracy. 


8 


III.  SCATTERING  OF  BENNETT-PROFILE  BEAMS 


A.  Preservation  of  Bennett  Profile 

On  the  basis  of  a  Fokker-Planck  treatment  of  scattering,  Lee4 
has  shown  that  the  Bennett  equilibrium  is  a  self-similar  solution  to  Eq. 
(2)  with  the  Bennett  radius  obeying  the  Nordsieck  equation,  Eq.  (6).  It  is 
desirable  to  know  if  this  important  result  carries  through  when  Moliere 
scattering  is  employed.  To  answer  this,  we  inserted  the  scattering  algo¬ 
rithm  described  in  Sec.  II  into  the  one-dimensional  particle  simulation 
code  SCATR.  SCATR  simulates  the  behavior  of  an  infinitely  long,  axisym- 
metric  self-pinched  beam  in  the  high-energy  paraxial  approximation  (v  =  c, 
v/y  <  1).  The  equations  which  SCATR  solves  are  (1)  the  Lorentz  equations 
for  the  transverse  motion  of  the  individual  particles  and  (2)  Ampere's  law 
for  the  self-magnetic  field. 

The  results  of  simulating  a  10  kA,  50  MeV  beam  with  an  initial 

Bennett  radius  of  0.5  cm  are  shown  in  Fig.  3.  In  Fig.  3(a)  the  mean  radius 

(the  radius  enclosing  half  of  the  total  current)  is  plotted,  and  is  seen  to 
e-fold  about  three  times  during  the  run.  The  fact  that  the  plotted  line  is 
straight  means  that  the  radius  is  growing  exponentially  with  a  constant 
Nordsieck  length  Lr.  In  the  Fokker-Planck  treatment  of  Eq.  (2),  this 
behavior  is  predicted  due  to  the  facts  that  (1)  energy  is  scattered  into 
the  transverse  direction  at  a  constant  rate  and  (2)  the  scattered  distribu¬ 
tion  is  Gaussian.  The  numerical  results  suggest  that  these  conditions  are 
at  least  approximately  fulfilled  when  Moliere  scattering  is  applied.  In 
Fig.  3(b)  we  compare  the  transverse  distribution  functions  of  the  beam, 
after  it  has  propagated  three  Nordsieck  lengths,  with  an  exact  Gaussian. 
The  agreement  is  good  except  in  the  tail.  Initially,  the  distribution 
functions  is  exactly  Gaussian,  as  required  for  a  Bennett  equilibrium,  and 
the  non-Gaussian  tail  Is  observed  to  grow  slowly.  At  the  stage  reached  in 

Fig.  3,  however,  it  appears  to  have  no  effect  on  the  spatial  current 


distribution,  shown  in  Fig.  3(c).  It  seems  reasonable  to  assume  that  the 
tail  will  have  little  practical  significance  for  beam  propagation:  it 
grows  up  so  slowly  that  the  beam  will  have  expanded  to  very  large  radius  by 
the  time  a  significant  number  of  particle  reside  in  the  tail.  As  the 
radius  increases,  the  confining  magnetic  field  weakens  and  the  state  of 
quasi-static  equilibrium  disappears,  allowing  the  particles  to  freely 
expand. 

The  presence  of  the  long  Rutherford  tail  on  the  scattering  dis¬ 
tribution  has  another  noteworthy  effect.  A  non-negligible  number  of  par¬ 
ticles  undergo  large  angle  scatterings  which  kick  them  clear  of  the  beam. 
If  a  particle  reaching  a  radial  position  greater  than  ten  Bennett  radii  is 
considered  lost  from  the  beam,  we  find  that  2-3%  of  the  beam  current  is 
lost  per  Nordsieck  length.  This  loss  rate  is,  of  course,  negligible  for 
practical  purposes. 

B.  Evaluation  of  Nordsieck  Lengths 

Equation  (7)  shows  that  for  a  Bennett  beam,  the  Nordsieck  length 
as  usually  defined  depends  only  on  the  beam  power  vy  (since  y2a2  is  a  con¬ 
stant).  Using  SCATR  we  have  found  this  to  be  the  case  with  Moliere  scat¬ 
tering  also.  The  Nordsieck  lengths  obtained  from  SCATR  are  compared  to 
those  from  Eq.  (7)  in  Fig.  4.  The  Rossi-Greisen  value  for  emax  was  used 
in  Eq.  (7).  The  implications  for  some  present-day  beam  machines  are  shown: 
for  PHERMEX,  the  SCATR  Nordsieck  length  is  twice  as  long  as  that  predicted 
by  Eq.  (7);  for  ATA,  the  SCATR  length  is  40%  longer  than  the  Eq.  (7) 
result.  When  vy  reaches  a  value  of  104,  there  is  little  difference  between 
the  SCATR  and  Rossi-Greisen  results.  Despite  appearances,  the  SCATR 
Nordsieck  length  is  not  exactly  linear  in  vy  (see  following  subsection). 


C.  Analytic  Expression  for  Nordsieck  Length 


The  envelope  equation  for  a  self-pinched  beam3  yields  a  general 
expression  for  the  Nordsieck  length. 


2vYm2c3 

<S>' 


(14) 


where  <dp^2/dt>  is  the  rate  at  which  the  mean  square  transverse  momentum  of 
the  beam  changes  due  to  scattering.  As  noted  in  subsection  A,  Moliere 
scattering  yields  (1)  an  approximately  constant  value  for  <dp^2/dt>  and  (2) 
an  approximately  Gaussian  distribution  in  transverse  momentum.  Thus,  to 
compute  the  value  of  <dp^2/dt>,  we  throw  away  the  non-Gaussian  corrections 
in  Mol i ere 's  Eq.  (11).  We  obtain. 


1 

T7 

me 


ny2a2  (vB  +  vt  $L)  (15) 


From  Eq.  (9)  we  obtain 

dB  1  B 

dt  ~  t  (B  -  1) 


(16) 


Substituting  Eqs.  (15)  and  (16)  into  Eq.  (14),  we  obtain 


2w  (B  -  1) 

,2  2.2 

ny  a  o 


(17) 


The  parameter  B  is  a  function  of  the  transverse  beam  temperature,  which 
tends  to  remain  constant  during  quasi-static  expansion  of  the  beam.**  This 
is  because  the  increase  In  transverse  energy  due  to  scattering  is  balanced 
by  the  work  done  In  expanding  against  the  magnetic  field.  To  compute  B  we 
use  the  expression  for  the  transverse  temperature  of  a  Bennett  beam,  i.e., 
vth2  ■  v/y ,  where  v*h  is  the  thermal  velocity,  so  that 
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(18) 


(19) 


(20) 


If  one  substitutes  values  for  a  and  n  appropriate  to  air  into  Eq.  (20)  one 
obtains 

lN  =  7*14  ^  meter- atm  .  (21) 

Equation  (21)  agrees  to  within  5%  with  the  results  of  SCATR  for  vy  >  50  and 
within  about  10%  for  vy  <  50,  as  shown  in  Table  1.  The  better  agreement  at 
larger  values  of  vy  is  to  be  expected  when  one  considers  that  corrections 
to  the  Gaussian  in  Moliere's  theory  enter  to  order  1/B. 
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TABLE  1.  Comparison  of  Nordsieck  lengths  obtained  from  SCATR  with  those 
obtained  from  Eqs.  (21)  and  (7).  The  value  of  the  Moliere 
expansion  parameter  B  used  in  Eq.  (21)  is  also  tabulated. 


Nordsieck  Length 

(meter-atm) 

VY 

SCATR 

Eq.  (21) 

Eq.  (7) 

B 

0.68  (PHERMEX) 

0.59 

0.55 

0.27 

6.6 

5.88  (ETA) 

3.9 

3.9 

2.3 

8.8 

8.0  (PHERMEX  UPGRADE) 

5.7 

5.1 

3.1 

9.0 

58.8  (ATA) 

33.4 

32 

23.2 

11.1 

588 

284 

272 

232 

13.4 

5880 

2404 

2367 

2320 

15.7 

13 


IV.  COMPARISONS  WITH  EXPERIMENTS 


A.  PHERMEX 

Our  work  on  small  angle  multiple  scattering  of  particle  beams 
originally  was  undertaken  in  support  of  a  series  of  propagation  experiments 
at  the  PHERMEX  electron  beam  facility  at  Los  Alamos  National  Laboratory. 
At  present,  PHERMEX  produces  a  200-300  A,  21  MeV  beam  with  a  1  nin  diameter 
beam.  This  beam  was  Injected  into  a  balloon  filled  with  helium  at  530 
torr,  and  the  spot  size  was  measured  as  a  function  of  propagation  dist¬ 
ance.  The  results  are  shown  in  Fig.  5.  The  limited  data  available  agrees 
rather  well  with  the  Moliere  scattering  result.  Additional  comparisons 
with  PHERMEX  will  be  made  early  next  year  at  an  order-of-magnitude  larger 
beam  power. 

B.  ASTR0N 

The  ASTR0N  beam  at  Lawrence  Livermore  National  Laboratory  was 
used  in  a  study  of  beam  expansion  due  to  scattering  some  years  ago.4*8 
In  one  experiment,  the  beam  was  injected  into  a  tank  of  gas  of  variable 
pressure.  The  beam  radius  was  measured  two  meters  downstream  as  a  function 
of  the  gas  pressure.  Because  the  ASTR0N  beam  has  a  low  power,  it 
Nordsleck-expands  only  for  a  short  while  and  then  goes  over  into  the  free 
expansion  regime.  As  a  result,  when  we  simulated  the  experiment  using 
SCATR,  we  had  to  use  a  "plural"  scattering  algorithm,  instead  of  a  multiple 
scattering  algorithm  as  discussed  in  Sec.  IIC.  The  experimental  results 
for  propagation  in  nitrogen  and  argon  are  compared  to  the  SCATR  results  in 
Fig.  6.  The  agreement  Is  good  for  both  gases,  and  that  for  argon  is  an 
Improvement  on  the  agreement  previously  reported  by  Lee.4  In  Refs.  4  and 
8,  the  value  of  emax  used  is  that  due  to  Williams9  and  is  approximately 
equal  to  (see  paragraph  following  Eq.  (11)).  However,  these  two  ref¬ 
erences  assume  that  emax  is  a  function  of  propagation  distance,  which  is 
true  for  noninteracting  particles,  but  not  for  a  quasi-static  self-pinched 
beam.  For  the  latter,  emax  *  (hta2B)1/2  is  a  constant,  with  B  given  by 


V.  EVOLUTION  OF  AN  ANNULAR  BEAM 


It  has  been  shown  by  Lee4  that  under  certain  restrictions  a 
self-pinched  beam  near  equilibrium  tends  to  assume  a  Bennett  profile  due  to 
small  angle  scattering.  The  rate  at  which  this  occurs  was  not  computed, 
however,  and  to  do  so  analytically  would  be  laborious.  SCATR  is  a  useful 
tool  for  addressing  this  problem.  We  chose  to  investigate  the  evolution  of 
a  beam  which  starts  out  far  from  being  a  Bennett  equilibrium  viz.,  a  rotat¬ 
ing  thin  annular  beam.  Like  the  Bennett  profile  beams  we  have  previously 
examined,  the  beam  is  assumed  to  be  charge-neutral ,  nondiamagnetic,  and 
self-pinched.  The  equilibrium  current,  magnetic  field,  and  transverse  vel¬ 
ocity  profiles  are  shown  in  Fig.  7.  The  value  of  vy  for  the  beam  is  50, 
which  gives  a  predicted  Nordsieck  length  of  27.5  m-atm  in  air  [from  Eq. 
(21)].  When  the  beam  propagates,  we  find  that  it  expands  slower  than  the 
predicted  rate  for  about  the  first  Nordsieck  length,  and  thereafter  expands 
at  the  predicted  rate,  as  shown  in  Fig.  8(a).  If  we  look  at  what  is  hap¬ 
pening  to  the  structure  of  the  beam  during  this  time,  we  find  that  its 
annular  shape  is  rapidly  destroyed.  Figure  8(b)  shows  that  in  less  than  a 
Nordsieck  length,  the  beam  has  completely  filled  in  and  its  transverse  vel¬ 
ocity  profile  is  close  to  being  Gaussian.  At  this  point,  only  the  rotation 
of  the  beam  is  preventing  the  current  profile  from  looking  more  like  a 
Bennett  profile.  As  the  beam  expands,  the  rotation  slows  down  due  to  con¬ 
servation  of  annular  momentum  and  the  beam  profile  becomes  very  close  to 
being  Bennett-like  (Fig.  8(c)).  On  the  basis  of  these  results,  we  think 
that  It  Is  reasonable  to  assume  that  any  self-pinched  high  energy  beam  with 
low  v/y  will  acquire  a  Bennett  profile  in  a  distance  on  the  order  of  one 
Nordsieck  length. 
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VI.  SUMMARY  AND  CONCLUSIONS 


We  have  carried  out  an  Investigation  of  multiple  small  angle 
scattering  of  particle  beams  based  on  the  Moliere  equations.  Our  results 
for  Bennett  equilibria  obtained  from  a  simple  particle  simulation  code 
SCATR  give  Nordsieck  lengths  significantly  longer  than  those  calculated 
from  commonly  used  expressions  based  on  the  less  accurate  Rossi-Greisen 
formalism.  On  the  basis  of  our  numerical  results,  we  are  able  to  derive  an 
improved  analytic  expression  for  the  Nordsieck  length  which  agrees  to 
within  5-10%  with  the  numerical  results,  depending  on  the  beam  power. 

We  obtain  good  agreement  with  experimental  data  on  beam  expan¬ 
sion  from  the  PHERMEX  and  ASTRON  electron  beams.  For  PHERMEX  we  predicted 
a  Nordsieck  length  which  was  twice  as  long  as  that  obtained  from  a  commonly 
used  expression.  The  experimental  evidence  supports  our  result.  For  the 
ASTRON  data,  the  agreement  is  significantly  better  than  that  previously 
obtained  by  other  workers. 

Finally,  we  have  studied  the  evolution  of  a  non-Bennett  equilib¬ 
rium  namely,  a  thin  annular  beam,  under  the  influence  of  small  angle  scat¬ 
tering.  In  a  distance  on  the  order  of  the  Nordsieck  length,  the  beam 
assumes  a  Bennett  profile  in  current  and  a  Gaussian  distribution  in  trans¬ 
verse  momentum.  Based  on  this  example,  we  tentatively  conclude  that  all 
self-pinched  beam  equilibria  with  v  »  c  and  v/y  «  1,  will  evolve  to  a 
Bennett  distribution  in  roughly  one  Nordsieck  length. 
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Figure  1.  Illustration  of  disparity  between  the  Moliere  and  Rossl-Greisen 
predictions  for  the  mean  square  transverse  momentum  of  a  beam  of 
noninteracting  particles  as  a  function  of  propagation  length. 
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Figure  4.  Nordsieck  length  as  a  function  of  beam  power.  Over  a  wide  range 
of  power  the  Nordsieck  length  obtained  from  SCATR  is  signifi¬ 
cantly  larger  than  that  obtained  from  Eq.  (7)  using  the  Rossi- 
Greisen  efflax.  The  implication  for  some  particular  electron 
beams  is  indicated. 


(b)  are  radial  profiles  of  the  current  density  and  magnetic 
field,  respectively,  while  part  (c)  shows  the  transverse  momen 
turn  distribution  function. 


